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We discuss progress in SM two-loop calculations of the pole position of the massive gauge-boson propagators. 



1. INTRODUCTION 

An important type of "universal" (process in- 
dependent) corrections relevant for complete elec- 
troweak 2-loop calculations of physical processes 
are the contributions from on-shell gauge boson 
self-energies which incorporate the relation be- 
tween bare, MS and on-shell (pole) masses (2-loop 
renormalization constant in on-shell scheme). 
Such calculations are important to scrutinize the- 
oretical uncertainties which might obscure e.g. 
the indirect Higgs mass bounds obtained by the 
LEP experiments. Most interesting are 2-^2 
fermion processes, which in future eventually may 
be investigated with much higher precision at 
TESLA in case the GigaZ [ |lj option would be 
realized. Here, we outline the complete 2-loop 
SM calculation of the pole position of the mas- 
sive gauge-boson propagators presented in [ ||, || . 

The position of the pole sp and the wave- 
function renormalization constant Zi of the prop- 
agator of a massive gauge-boson are defined by 
the relation 
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der by order. Up to two loops the solution reads 

m 2 + n (i) + nP) + nWnW', (2) 
i-nW-nW'-nWnW", (3) 

where Il( L ) is the bare or MS -renormalized L- 
loop contribution to II, the prime (double prime) 
denotes the derivative (second derivative) with re- 
spect to p 2 at p 2 = m 2 . One of the remarkable 
properties of (0) is that the pole position is well- 
defined in terms of self-energy diagrams and its 
derivatives at momentum (square) equal to the 
bare or the MS mass (square) which, by con- 
struction, are real parameters. 

For the Z-boson the equation for the position 
of the pole is modified due to 7 — Z mixing 



s P 



m 2 z 



n~zz(sp) 



s P 



n 



77 (s P ) 







The 2-loop wave-function renormalization con- 
stant in this case is equal to 
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for s ~ sp where Tl(p 2 ,---) is the transversal 
part of the one-particle irreducible self-energy and 
TOq is the bare mass and we have adopted the 
standard parameterization of the pole spy = 
My — iMyTy, where Ty is the width of the un- 
stable gauge-boson. 

This equation can be solved perturbatively, or- 
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We would like to stress, that in order to obtain 
gauge invariant results the tadpole contributions 
have been included [ || . 

2. CALCULATIONS 

At the 2-loop level we have about 1000 1PI di- 
agrams each for the Z- and the W^-propagator. 
Our calculation is largely automatized: we use 
QGRAF [ H to generate the diagrams and then 
the C-program DIANA [ |(| to produce for each 
diagram the input suitable for our FORM pack- 
ages. 



2 



Two-loop propagator type diagrams with sev- 
eral masses can be reduced to a restricted set of 
so-called master-integrals by using a complete set 
of recurrence relations given in [ ffl] . We have used 
this approach only for the calculation of the 2- 
loop massless fermion corrections, where analyti- 
cal results for the master-integrals are available [ 
||. As compared to the existing calculation of 
the massless fermion contribution, performed in [ 
we apply Tarasov's recurrence relations [ ^ 
which allows us to reduce the number of master- 
integrals to a minimal set. The latter includes 
integrals which may be evaluated by using the 
package ON-SHELL2 [ [kJ plus the following- 
new prototypes: (using standard notation [00]) 

• ZZ: 

Fwowoo 7 Fqooow , Vhooz, Vwoow, Joo{h,w} ■ 

• WW: 

Fzowoo>FqoooZ) Vwooz, V{h,z}oow> Jqo{h,z} ■ 

We have worked out an independent analytical 
calculation of these master integrals by using a 
method developed in [ ||. One of the crucial 
points of our approach is the observation [ [llj 
that the analytical results for diagrams with two- 
massive cuts have a very simple form if written 
in terms of new variables. For diagrams with two 
equal internal masses m 2 and arbitrary external 

momentum p 2 the new variable is y = - — 



where z = p 2 /m 2 . When the external momen- 
tum belongs to one of the internal masses p 2 = 

WOi where 
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2 the new variable 



is X 



l + y/l-ix ' 

x = m 2 /M 2 . To speed up the numerical eval- 
uation of on-shell gauge-boson self-energies, an 
expansion in sin 2 Q\y may be used, which con- 
verges well and simplifies the results considerably. 
However, it is not a naive Taylor expansion due 
to the presence of threshold singularities of the 
form In 7 sin 2 9w ■ We therefore briefly outline the 
corresponding expansions for the relevant master- 
integrals by giving two examples. The integral 
Fwowoo [z = 1/ cos 2 6w] depends on the small 
parameter s = 1 — 1/z via the variable (see [EJ) 
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which we expand first and then substitute it into 
the analytical result in order to get the series ex- 
pansion for Fwowoo- The integral Fzowoo[x = 
cos 2 9w] depends on s via the variable 



X = 




13=0 



exp i 



2?r 



'71 ^ 

V" 3 n=0 



2n 



1 



and we proceed as in the previous case. How- 
ever, in all other cases, the master-integrals which 
show up are not know analytically so far. They 
can be calculated numerically, by using one of the 
approaches and/or programs developed for the 2- 
loop self-energies [ ^2|. The full list of master- 
integrals occurring in the 2-loop calculation of 
the pole-mass of the gauge-bosons (for the photon 
propagator see [ 13 ) reads: 



• ZZ: 

FwWWWH i FuttH , FzWHWZ i FzHHZZ i FzZHHH , 

Fwwwwz , Fztntt, Fwowot, FwtwtO: F t omw, 
VwZWW > VwHWW i Vh WWZ , VjffttZ , Vzwwh , 
VzttH , VtHtt , Vwotw 7 VtOWt 7 Vtztt , Vhhzz 7 Vzzzh , 
Jzww, Jhww, Jowwi Jotti Jztti Jhu, Jqwz-i 

JzHH, JoWt ■ 



• WW: 

Foottz, Fy/oztt, Fwomt , Fwwzzh , F\ 

FwWHHH , Fy/ZZWW , Fy/ZHWW 7 Fy\ 



^WWHZZ , 



vv ujj (-r 1 — vv vv *j ±± ■) — vy v\ li t-t z-i i 

■r WWHHH , ^WZZWW , Fy/ZHWW 7 Fy/HH WW 7 

Fwtzoo, Voowt, Voztt, VoHtti Vwhhh,Vwwwz, 
Vwwwh , Vwttz 7 VwttH , VwZZH 1 Vy/ZHZ , 
Vzotw 7 Vnotw 7 Vzwhw , Vhwzw , Vzwzw , 

VhWHW, JoWZ, JoWt, JqWH, JoZt, JoHt, JwHH, 
JwZZ, Jwtti JwZH ■ 

To keep control of gauge invariance we adopt the 
i?£ gauge with three different gauge parameters 
£,Wi £z and 67 ■ in niost cases exact analytic re- 
sults in terms of known functions are not avail- 
able. Thus, instead of working with the exact for- 
mulae (which only can be evaluated numerically, 
at present) we resort to some approximations, 
namely, we perform appropriate series expansions 
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in (small) mass ratios [ EBJ. For diagrams with 
several different masses it is possible that several 
small parameters are available. In this case we 
apply different asymptotic expansions (see [ llJJ) 
one after the other. Specifically we expand in the 
gauge parameters about £j = 1, in sin 2 8w and, 
for diagrams with Higgs or/and top-quark lines, 
in m v /m 2 H or/and my/m 2 . Numerical results 
are obtained using the packages ON-SHELL2 [ 
and TLAMM [ |IJ) (see [ § §). 

Renormalizing the pole-mass at the 2-loop 
level requires to calculate the 1-loop renormal- 
ization constants for all physical parameters 
(charge and masses), as well as the 2-loop mass- 
renormalization constant itself. Not needed are 
the wave-function and ghost sector renormaliza- 
tions. 

The full 2-loop relation between pole and MS 
parameters can be written in the form 



spy 



L v,o 



n 



(i) 



n 



(2) 
V.O 



n «n (1)/ 

1L v.o 11 v.o 



E 



(<Ko) (1) 



dm 2 



de 



(4) 



where the sum runs over all species of particles, 
(<5eo)^ and {5m 2 are the 1-loop counter- 
terms for the charge and physical masses in the 
MS -scheme. The derivatives in (||) correspond to 
the subtraction of sub-divergencies. The genuine 
2-loop mass counter-term is obtained by expand- 
ing Tfiy q in terms of the renormalized mass. 

Summing all 2-loop contributions we restore 
gauge invariance of the position of the com- 
plex pole before UV renormalization. After UV- 
renormalization the propagator pole is repre- 
sented in terms of finite amplitudes. For explicitc 
results we refer to [ ||, || . For the numerical eval- 
uation we were using the first six coefficients in 
the weak angle sin 2 6w and in the mass ratios 

2/2 
m V I m heavy 

3. RENORMALIZATION GROUP EQUA- 
TIONS 

We now consider the SM renormalization group 
(RG) equations. We will denote on-shell masses 
by capital M and MS masses by lowercase m. 
We adopt the following definitions for the RG 



functions: for all dimensionless coupling con- 
stants, like g, g' , g s , e, A, the /3-function is given 
by A* 2 5 — Pg an d for all mass parameters (a 
mass or the Higgs v.e.v. v) the anomalous dimen- 
sion 7 m 2 is given by fi 2 -J^ In m 2 = 7 m 2 . Using the 

fact that sp is RG- invariant: 1^-^iSp = 0, we are 
able to calculate the anomalous dimension of the 
gauge bosons masses from our finite results. At 
the same time, a typical relation between bare- 
and MS -masses has the form 



4,o = m v(p) (1 + 



k=l 

such that the RG functions may be calcu- 
lated directly from the UV counter-terms jy = 
Ei \9i~£^ Z V ' C? = 9,9s)- In addition, the 
UV counter-terms satisfy relations connecting the 
higher order poles with the lower order ones: 
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In the SM it is interesting to compare the RG 
equations calculated in broken phase with the 
ones obtained in the unbroken phase. Let us re- 
mind that at the tree-level the vacuum expecta- 
tion value v is given by v 2 = where m 2 and 
A are the parameters of the symmetric scalar po- 
tential. The masses of the gauge-bosons are equal 
to m 2 z = j(g 2 + g' 2 )v 2 and m 2 ^ = \g 2 v 2 , respec- 
tively. The fact that these relations are RG in- 
variant on the level of the bare quantities implies 
the relations 
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where the 2-loop RG functions /3 g , @ g > , (3\ , 7 m 2 
have been calculated in the unbroken phase in [ 
p"7| . We have verified in the MS scheme, that 
these relations are valid up to 2-loop order in the 
broken phase with the same RG functions. Thus 
the RG equations for the MS masses in the broken 
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theory can be written as 
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m w (n ) 
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rn 2 H ( M 2 ) - 2mV) , 
1^ 
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2 A(^ 2 , 

where y t is the top-quark Yukawa coupling (the 
other Yukawa couplings are kept zero). The MS 
Fermi constant 



G F (^) EE 



(9) 



(10) 



8m 2 v (ii 2 )sin 2 d w {^ 2 ) 
satisfies the following RG equation 

2 ® i A / 2\ #A 

The knowledge of the anomalous dimensions 
7v allow us to write expression for the pole po- 
sitions spy with explicit factorization of the RG 
logarithms. What we get is 
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where fi 2 In m 2 
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and 9 = 7 { / . 1) + 2(3 { g ] 
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In 



(1) (i i) 

contrast to QCD, j v and C v have non- 
polynomial structure in the massless coupling 
constants which originated from the tadpole con- 
tributions. Xy t x has been calculated long time 
ago [ [|] and Xy.2 are our results which, for the 
numerical evaluation, we approximated by finite 
series. The functions Xy.j we have represented 
in terms of MS parameters. We note that the 
amplitudes Xyj entering ([ll]) have no explicit ji 
dependence. 



4. SCHEME DEPENDENCE 

Our results reveal terms of unexpectedly high 
powers of the Higgs and the top-quark masses in 
(pi]), arising from 2-loop corrections. In fact the 
purely bosonic diagrams yield m A H /m v terms and 
the (t, b) quark-doublet (m t 3> mb) contributes 
rn^/[m H rny) power corrections. At a first glance, 
such terms contradict Veltman's screening theo- 
rem [ [l!| which states, that the L-loop Higgs de- 
pendence of a physical observable is bounded by 
(ttVj I ) l ~ 1 In m\ for large Higgs masses. How- 
ever, this theorem only applies to physical observ- 
ables like cross sections and asymmetries. If we 
consider quantities like Ar (which is an observ- 
able in the on-shell scheme) in the MS scheme the 
screening theorem does not hold in general. To 
illustrate this, let us compare the l-loop EW cor- 
rections to the Fermi constant [ Q in the on-shell 
and the MS schemca: 

G F = —— — — ll + Ar. 



\/2A/ 2 sin 2 6 
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so that 
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The l-loop correction Ar^ in the limit of a 
heavy Higgs boson has a logarithmic Higgs mass 
dependence only . In contrast, in the MS scheme, 
higher powers of the Higgs mass are showing up, 



because 



Mr, 



1 contributes the term 



Ar 
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MS 
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while all other corrections exhibit logarithmic be- 
havior in the Higgs mass, onlyB. Our relation (|ll"| ) 
between parameters of two different schemes does 
not relate physical observable but it is very im- 
portant for analysis of the uncertainties coming 



2 We should mention that Ar^g introduced in (|L4|) is dif- 
ferent from Ar defined via, a hybrid scheme (couplings MS 
, masses on-shell) in [ bo|] . 

3 The hybrid quantity Ar^ 1 ' does not show extra powers 
of the Higgs mass since masses are kept on-shell. 
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from higher order effects [Um. In particular, our 
estimations of the 2-loop boson contribution to 
Ar [ H is very close to the results of real calcula- 
tions [||, f§. 

5. CONCLUSION 

By an explicite 2-loop calculation we have 
shown that (to this order): 1. The position of 
the complex pole s p of a the gauge-boson (Z, W) 
propagator is a gauge invariant quantity after in- 
clusion of the Higgs tadpole contributions (see 
also [ [24| )• 2. The renormalized on-shell self- 
energies are infrared finite. This derives from 
the fact that within dimensional regularization, 
which allows to regularize UV and IR singulari- 
ties by the same e (e = (4— d)/2 — > 0) parameter, 
the singular 1 je terms are absent after UV renor- 
malization. 3. The inclusion of the tadpoles is 
important for the renormalization group invari- 
ance and for the gauge invariance of the param- 
eter renormalizations. 4. By our calculation we 
have proven that the MS renormalization scheme 
is self consistent and works properly in case of un- 
stable particles. 5. Our results for the 2-loop mass 
renormalization constants in the on-shell and the 
MS scheme can be applied in calculations of phys- 
ical quantities in both of these schemes at the 2- 
loop level. Examples of such calculations, where 
results of our paper [ |) have been used, are the 
computations of the bosonic 2-loop contributions 
to the Muon life-time presented in [ . 
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